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Abstract. Parton-shower event generators that go beyond the collinear-ordering approximation 
at small x have so far included only gluon and valence quark channels at transverse momentum 
dependent level. We describe results of recent work to include effects of the sea-quark distribution 
with explicit dependence on the transverse quark-momentum. This sea-quark density is then applied 
to the description of forward Z -production. The qq* — > Z matrix element (with one off-shell quark) 
is calculated in an explicit gauge invariant way, making use of high energy factorization. The kj- 
factorized result has been implemented into the CCFM Monte-Carlo CASCADE and a numerical 
comparison with the qg* — > Zq matrix element has been carried out. 

PACS: 12.38.-t,12.38.Cy,12.39.St, 14.70.Hp 

INTRODUCTION 

Transverse momentum dependent parton distribution arise at small x naturally as a con- 
sequence of high energy factorization and BFKL-evolution [1]. A formulation of high 
energy factorization which is in accordance with conventional collinear factorization is 
provided by kj -factorization as defined in [2, 3]: at high center of mass energies, the 
resummation of high energy logarithms (BFKL) can be brought into a form consistent 
with conventional collinear resummation (DGLAP). The CCFM evolution equation [4], 
which interpolates between DGLAP and BFKL evolution, allows then for a Monte Carlo 
implementation of kj -factorization, which is realized by the Monte Carlo event genera- 
tor CASCADE [5]. 

Based on the principle of color coherence, the CCFM parton shower describes only 
the emission of gluons, while real quarks emissions are left aside. Note that such an 
approximation is in principle justified as enhanced regions of phase space for x — > 
1 and x — > are dominated by gluonic dynamics at leading logarithmic level. For 
the description of unintegrated parton densities this implies that the CCFM evolution 
describes only the distinct evolution of unintegrated gluon and valence quarks [6] while 
non-diagonal transitions between quarks and gluons are absent. On the other hand, it 
is necessary to include quark emissions into the parton shower to take into account 
subleading effects. Also the determination of kj -dependent hard matrix elements is 
affected by the absence of quark emissions: gluon-to-quark splittings needed for the 
generation of seaquarks are to be included into matrix elements, which complicates 
their structure, see e.g. [7]. In order to supplement the CCFM evolution with quarks 



we present in this contribution a definition of a partonic matrix element involving off- 
shell quarks and the definition of an off-shell sea-quark density. For the latter we restrict 
here to the case where the gluon-to-quark splitting occurs as the last evolution step. To 
test our procedure numerically, we apply it to the production of a Z-boson in the forward 
direction at LHC energies. For a detailed discussion we refer to [8, 9]. 



DEFINITION OF A q T -DEPENDENT SEA-QUARK 

Let us take the approach of generalizing kj -factorization to the case of quarks by mim- 
icking the already existing gluonic case. There off-shell gauge invariance is ensured 
through a reformulation of QCD at high center of mass energies in terms of effec- 
tive degrees of freedom, reggeized gluons. The latter coincide in their on- shell limit 
with conventional collinear QCD gluons, while one uses for the general off-shell case 
effective vertices which contain additional induced terms ensuring off-shell gauge in- 
variance 1 . In complete analogy one can construct a reggeized quark formalism for the 
description of the high energy limit of scattering amplitudes with quark exchange in the 
/-channel [12]. As (reggeized) quark exchanges are suppressed by powers of s in com- 
parison to (reggeized) gluon exchanges, they generally do not occur in the high energy 
resummation of total cross-sections. They can however be used as a starting point for the 
construction of an off-shell factorization of matrix elements which are limited to quark 
exchange in the /-channel. This is the case for the g*q — >■ Zq matrix element for which 
we construct [8] an off-shell factorization into an off-shell matrix element q*q — > Z and 
a corresponding gluon-to-quark splitting function. Making use of the reggeized quark 
formalism, we arrive at the off-shell partonic cross-section qq* — > Z, 

a{x l x 2 sMW) = V2G F M 2 {V 2 +A 2 q )^8{x l x 2 s-q 2 -M 2 ). (1) 

Here x\ y/s and x 2 \/s are the light-cone momenta of the two incoming quarks, while q 
denotes the transverse momentum of the off-shell quark; Mz denotes the mass of the 
Z-boson, while \flGpM\(y 2 + A 2 ) describes the coupling of the Z-boson to the quark. 
Eq. (1) agrees for q 2 — > with the qq^rZ cross-section, while the splitting function 
is obtained as a pure color factor if the reggeized quark formalism is taken literally 
— a crude approximation which can be traced back to non-conservation of energy at 
the level of high-energy leading logarithms. It is possible however [8] to supply energy 
conservation while maintaining all desired gauge invariance properties. Doing so we 
arrive at the kj -dependent [3] gluon-to-quark splitting function 
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where k and q denote transverse momenta of the off- shell gluon and quark respectively, 
while z denotes the fraction of gluon light cone momentum which is carried on by the 



For an effective action approach to the reggeized gluon formalism see [10, 1 1]. 



/-channel quark. In the on-shell limit k 2 — > eq. (2) reduces to the collinear DGLAP 
splitting function Tr[z 2 + (1 — z) 2 } at lowest order; in the high-energy limit however it 
factorizes correctly the finite k dependence to all orders. At the next-to-leading logarith- 
mic accuracy a s (a s hxx) n , the qj -dependent sea-quark density is defined by 




with k 2 ^ = \l 2 jz — q 2 /(z(l — zj). For the scale \x of the unintegrated gluon den- 
sity, we will investigate two possible choices. One choice, closer in spirit to inclu- 
sive calculations [3], is to set simply p, 2 = /i 2 . A choice based on angular ordering, 
which is more natural from the point of view of the CCFM evolution, is given by 
fi 2 = q 2 /(l — z) 2 + k 2 /(l — z). The complete gy-factorized cross-section for forward 
Z-production is given by 

°%~4z\s,Ml) =£ fdx, fdq 2 f dx 2 a! iq ^ z ^(x u q 2 ^ 2 )Qf(x 2 ^ 2 ) (4) 
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with Qj (x2, li 2 ) the valence quark distribution, implemented as in the first paper of [6]. 



NUMERICAL TEST OF THE OFF-SHELL FACTORIZATION 
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FIGURE 1. 

We next compare numerically the factorized expression eq. (4) with the result of 
implementing the g*q — > Zq off-shell cross-section [8, 13] in the CASCADE Monte Carlo 
generator [5]. The results are summarized in fig. 1 and fig. 2. For the factorized result 
both choices for the hard scale p, 2 have been implemented. The differences between the 
matrix element (ME) and the factorized results are discussed in [9]. We note that the 
factorized result contains the /-channel contribution, while the s-channel contribution 
is sub-leading both in the collinear and high energy limit. The factorized expression 
also contains a kinematical approximation [8, 12] in comparison to the full g*q — > Zq 
cross-section, which leads to a small, but finite correction. 

In summary, we achieved an off-shell factorization of the g*q — > Zq process, where the 
factorized expression can both analytically and numerically be shown to coincide with 
the full g*q — > Zq cross-section up to a finite remainder, sub-leading in the collinear and 
the high energy limit. 
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FIGURE 2. 
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